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Introduction 

This paper begins a series on the estimation within satisfactory tolerances of the 
number of character weighted lattice points in natural curved regions of Euclidean 
spaces and applications of the results and methods to the study of zeta- and L-functions. 
The study of the lattice points in a circle presents a model for the subsequent treatment 
of the Dirichlet and other divisor problems. 

C.F. Gauss [Gausl] [Gaus2] showed that the the number of lattice points in circles is 
the area, irt , up to a tolerance of order at most the circumference, i.e. the square root, 
0{^/i) . Sierpinski in 1904, following work of Voronoi [Vor] on the Dirichlet Divisor 
Problem, lowered the exponent in Gauss's result to 1/3 [Sierpl] [Sierp2]. Further 
improvements that have been obtained will be described below. This paper provides a 
proof of the conjectured optimal error bound, 0(t4+'^) . 

The present approach begins with a new asymptotic expression (9) for the number 
of lattice points in a circle. To obtain such an expression, one begins by embedding the 
integer lattice within a lattice refined by a factor Q that can be adjusted depending on 
the size of the circle. A convex polygon with vertices in the refined lattice approximates 
the circle. This number of lattice points can then be rewritten as a sum of standard 
trigonometric functions over the points of the refined lattice that lie within this poly- 
gon. An elementary generalized Euler-Maclaurin formula (the shortest possible such 
formula) for lattice sums in lattice polygons then provides an estimate of this sum in 
terms of certain integrals over the polygon. The lattice polytope structure permits such 
a formula in terms of Lebesgue measure. The integrals can in turn be approximated 
by integrals over circular regions; classically, such integrals yield Bessel functions. This 
procedure thus leads to an expression for the number of lattice points in a circle in 
terms of sums of Bessel functions or exponential sums. The expressions considered by 
Hardy, Landau and others have a similar form, but (9) differs in important features; 
for example, the sum is taken over lattice points in a square, whereas classically the 
summation pattern in asymptotic formulae for the lattice points in a circle retains the 
circular symmetry. 

The classical approach to any such sum next calls for application of Poisson sum- 
mation, followed by lots of estimates, with the goal of arriving at a situation in which 
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the method of stationary phase apphes to transform the original problem into a new 
one. It turns out that only a portion of the summation has relevance to proving the 
conjectured error estimate; the rest contributes less than the desired error. The present 
approach considers an expression S{t), defined in (10), that in fact contains a form of 
the crucial terms. The nature of the integrals to be summed over (m, n) suggests di- 
viding the domain of summation into regions near and far from the circle + = t . 
The method of (9) and what preceeds it plays a key role in the analysis of the sum over 
points near the circle, for example in the proof of Prop. 3.1. The final result, (17), 
is a complete analytic expression for S{t) for a given division of the domain into near 
and far regions. Upon comparison of different choices of this division depending on a 
small enough parameter rj > , some terms drop out and what remains provides the 
conjectured estimate in the error term in the circle problem. 

The history of our efforts related to these problems begins in the early '90's with 
results on general Euler-Maclaurin expansions, including with remainder, for lattice 
convex polytopes, [CS4 ] [CS2 ] [S], and their relations to certain reciprocity laws for 
some generalized Dedekind sums, extending known connections between lattice prob- 
lems and certain Dedekind sums [P]. The possibility of such applications provided an 
impetus for investigations of Eulcr-Maclaurin expansions and possible number theory 
connections. Our approach generalized methods we had earlier used to compute all co- 
efficients of the Hilbert-Ehrhardt polynomials that count the number of lattice points 
of convex lattice polytopes [CSl]. These methods relied on our computations of the 
Todd classes of toric varieties, in turn based on our general formulae to compute (many 
kinds of) characteristic classes of singular varieties in various topological and algebraic 
settings [CSl ] [CS3 ] (see also [CMS1,2], [CLMS]). Higher dimensional analogues of 
the Euler-MacLaurin formula will appear in subsequent papers in this series. Some 
other valuable modern contributions on lattice points and Euler-MacLaurin formulae 
from differing perspectives include [BV] [DR] [K] [Kh] [KK] [KSW] [M] . 

A summary of earlier work: A shorter proof of the Sierpinski estimate was given by 
Landau [Lan2] [LanS] in 1912, after making rigorous a heuristic geometrical method 
of Pfeiffer [Pfe] of 1886. Hardy [H1915] [H1916a] and Landau showed in 1915 that the 
error bound could not be lowered to 0(1^/"^) . The Hardy-Landau formula [L] for lattice 
points in circular regions constituted the basis of subsequent efforts. New methods 
of estimating exponential sums were introduced by Hermann Weyl, and by Van der 
Corput [VdCl] [VdC2], who in 1923 gave exponents slightly below 1/3 for Gauss's circle 
problem. His estimate of 33/100 was subsequently improved in work of his student L.W. 
Nieland [Nie] to 27/82= 0.329268... ,which Van der Corput himseff had obtained for the 
Dirichlet Divisor Problem. Littlewood and Walfisz [LW] shortly thereafter lowered the 
estimate, initially obtaining 37/112 = 0.33035..., followed by further small improvements 
[Wal] [Wa2] [Wa3]. Further improvements were obtained by Vinogradov 1932 to 34/106 
= 0.32075... , and by Titchmarsh in 1934 [Ti]. Other results obtained include Hua 
Loo Keng [Kengl] in 1942 of 13/40 = 0.325; Yin Wen Lin [Lin] in 1962 and Chen 
Jing Run [Run] in 1963 of 12/37 = 0.324324...; and Nowak [Now] in 1984 of 35/108 
= 0.324074... . The theory of exponential pairs gave a systematic account of some 
approaches to exponential sums; see [GK]. G. Kolesnik [Kol] obtained the exponent 
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bound of 139/429 = 0.324009324... .Using methods introduced by E. Bombieri and 
H. Iwaniec [BI] and further development, Iwaniec and C.J.Mazzochi [IM] obtained in 
1988 the then best bound of 7/22 = 0.3181818... M. Huxley [Hul990] improved this 
estimate to 46/146 = .31507..., and then slightly further in [Hul996] [Hu2001]. 

The present paper represents a considerable simplification and shortening of a ver- 
sion of a several years ago. We benefitted from comments from colleagues at various 
points, including Enrico Bombieri for surveying for us the theory of exponential pairs, 
Peter Sarnak on several occasions, David Kazhdan, M. Huxley, V. Guillemin and S. 
Sternberg, Steve Shatz, and Ed Miller for a careful reading of this version in a slighly 
longer form. 

The research is partially by a grant from DARPA. 

1. An Approximate Euler MacLaurin Formula in the Plane 

Let 7? C be the lattice of integer points in the plane. Let 

I t e Z 

Let TTi , i = 1, 2 be projection on the first and second coordinates. If v = (f i, ^2) , let 
V* = (f 2, and = {—V2, vi) . Let S be the sum over all integers of Dirac 5-functions 
concentrated at integer points. Specifically, it will be convenient, as pointed out by S. 
Sternberg, to use a Gaussian approximation, i.e. to take 

S — lim 5e , 




If = 1 — 5e and V'i,e(0) = 0, then ipi{t) = limg^o+ V'i,e • Therefore, we use the 
notation 

For E a line segment in the plane with integral cndpoints, we assign normalized 
Lebesgue measure so that the interval between two adjacent lattice points has length 
one. If the endpoints are (a, c) and {b, d) , let A be the greatest common divisor of c — a 
and h-d and let (mi, 777,2) = ^) • Let 

Lix) = {x — a) + c. 
mi 

Then integration over E with respect to normalized Lesbesgue measure satisfies 

// = — f f{x,L{x))dx=— f f{L-\y),y)dy. 
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There is no assumption needed on the order relation of a and 6 or c and d . The first 
equahty assumes the segment is not vertical, the second that it is not horizontal. This 
holds even for / a bounded function on L with finitely many discontinuities at each 
of which it posesses (finite) one sided limits /"*" and /~, and if such an / is the limit 
/ = limg^o+ fe of uniformly bounded continous functions, 



/ / = hm / /, . 

Je e^0+ Je 



On the other hand, let 



Note that 



1 

«L = - lim /'5eWl,e(^t) 

1 1 

Jo Jo 



2' 



as {ipi^eit)] < I • Hence aL <\ - Then for / continuous we have 

lim / /• (^eOTTi) • (V'i,e0 7r2) = lim — /" f{x,Lx)5^{x)il:i^^{Lx) 
e^0+ Je e-*o+ mi J„ 



1 

mi 



^ /(m, Lm)'ip{Lm) + aL{f{b, Lb) - f{a, La)) 

me(a,6)nZ 



Briefly, this can be seen as follows. Let 77 > be small. Then 



lim / f(Lx)5e{x)'ipi^e{Lx) = lim 



0+ 



a+rj 



f{Lx)6e{x)'lpi^e{Lx) 



/m+T] f-b 
f{Lx)5,{x)iPi,,{Lx)+ / f{Lx)5,{x)iPi,,{Lx) 

IfmeZ,a<m<6 and Lm is not an integer, then 

fm+r] 



lim 



f{x, Lx)6e{x)il;i^e{Lx) = f{m, Lm)ipi(Lx) 



m — rj 



On the other hand, 5e{—t) = 5e{t) and ipi^ei—t) = —ipe,i{'t) ■ Hence, if Lm e Z , 



/ f{x,Lx)6,{x)i;i,,{Lx)^ {f{m+t, L{m+t))-f{m-t, L{m-t)))6,{t)i;i,,{^J)dt. 

Jm—r) Jo 

Hence, given rj' > , there is an rjm such that 
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{f{x+t,L{x+t))-f{x-t,L{x-tm^,,{^^t)dt < Tj' 5,{t) < V 



Similarly, for example, at the left endpoint 

r^'' fix,Lx)S,ix)il;i,,iLx) = [\fia + t,L{a + t))S,it)il;^,,i^J)dt 

Ja Jo 

= f\f{a,La)5S)M^t)dt+ f\f{a + t, L{a + t) - f{a, La))5,{t)iP,,,{^t)dt , 
Jo Jo 

The first expression on the right converges to —aLf{a, La) , and by the same argument 
as for the interior integer points, there is 77 = ?7a > such the second expression will 
be less than r]' in size. Similary for the right endpoint. It follows that 

lim / / ■ (5e o TTi) ■ (V'l.e o 7:2) = hm — / f{x, Lx)6e{x)ipi^^{Lx) 

= —\ E /(™' Lm)i;{Lm) + aL{f{b, Lb) - f{a, La)) + E 
me[a,6]nz 

\E\ < {h — a + l)r]' , for any rj' > . 
We will use the notation 

/ / • ((^O TTi) • (V'l O 7r2) = lim / f ■ (Se O m) ■ (ipi^e ° 7^2) ■ 

Je ^^0+ Je 

(Warning: The value of this expression depends on the specific choice of approximating 
continuous functions.) 

Similarly, with the variables interchanged, 

/ /(V'l O TTl) • (S o 772) 

Je 

— \ Yl n)HL-'ri) + a^-i (/(L-M, d) - /(L-^c, c))l . 

1712 L J 



ne[c,d]nZ 



Integration in the plane will be with respect to Lebesgue measure with the unit 
square having unit area. Using = 1 — S , we also write 



L>(n,m) ((V'l o ^1) ■ (V'l ° ^2)) = n(l - 5 o TTl) • {ipi o 7^2) + m(^l o TTl) • (1 - 5 O 7T2) 
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(1.1) Theorem. Let A be a convex polytope in with vertex points in . Let 
Ei,...,Ek be the edges of A , listed going counter-clockwise around the boundary, say, 
and let nj e be the minimal integral outward vector normal to Ej . Let Ej have 
normalized Lebesgue measure so that the interval between two adjacent lattice point has 
measure one. Let f be a function on . Let 

T(A,/)= [ f + [ /.•(V'lOTri) + / /.•(V'lOTra) 

J A J A J A 

+ / /xj/ • (V'l °7ri)(V'i o7r2) 

J A 

1 ^ f 

i=i ''^i 



Then 



fix,y) - T(A,/) < Yl 



{x,y)£Z^nA 



(x,y)£Z^ndA 



To prove this result, let 



: Z^ n A ^ (0, 1] 



be a weighting of the lattice points in the polytope defined as follows: If (m, n) is in 
the interior of A , then w{m, n) = 1 ; if (m, n) lies on an edge but is not a vertex, then 
w{m, n) = ^ ; and if (m, n) is a vertex, then w{m, n) is the proportion of the area of 
small circle centered at (m, n) that lies within A . Then 

Y ^{^^y)f{^^y) ^ f ■ {S oni) ■ {5 0712) . 

N^r»2^ A J A 



{x,y)ez^nA 



We may write 



/ / •a07ri)-(5,07r2)= / / - / f-{l-5,OTT^)-! /-(l-^eOTTs) 
J A J A J A J A 

+ / /•(l-5,o7ri)-(l-5,o7r2). 

We first integrate by parts with respect to y . Let For Ej a non- vertical edge, let 
{aj,Cj) and {bj,dj) be the endpoints of Ej , with aj < bj,. Let \j be the greatest 
common divisor of bj — aj and dj — Cj , let nij = (mj^i, 777.^^2) = ( , '^^ ) and let 
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L^lx) = — —ix — a,- 1 + Cj 

■' TOj 1 



By integration by parts, 



/ /■(l-5e07r2) = - / fy{i;i,,on2)+^^^ f'^ f{Lj{x))i;,,,{L,{x))dx . 

J A J A J TTinij J a. 

The sum is over no n- vertical edges. The sign keeps track of whether mj- and rij are 
pointing in the same or opposite directions. Hence 

k 

/ / • (1 - 5e O 7r2) = - / fy {'lpl,e O 7r2) + ^ n2nj / / • (•(/'i^e O 7^2) ■ 
J A J A ^-^j JEj 

Similarly 

k 

/ / • (1 - 5e OTTi) = - / /a: ■ (V'l.e O TTi) + ^TTiIlj / / " ('i/'l.e O TTi) . 
J A J A JEj 

Next, mutiply these equations by 1 — (5e o tti and 1 — d^o n2 and integrate the area 
integrals with respect to x and y , respectively. The average of these two results is 

/ /•(l-5e07ri)-(l-5e07r2) = 
J A 

k n k ft 

-^niXlj / • (V'l.e OTTi) • (1 - 5e 7r2) - ^7r2nj / / • (1 - 5e OTTi) • (V'l.e 7r2) 



The last two terms can be combined as 



^TTiUj /• (V'l.eOTTi) • (1 -5e0 7r2) + 5^7r2nj- / / ' (1 - O TTi ) • (V'l.e O 772) 

= / / ■ ^n* ((V'l.e O TTi) • (V'l.e O 772)) . 

and then the product rule for derivatives and letting e go to zero finish the proof. 

Notes: 1. Convexity is not being used in any essential way. 
2. Clearly there is an exact formula 
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wix,y)fix,y) = T{A,f) 

(a;,y)eZ2nA 

for the weighted sum introduced in the proof. This is the first, and shortest, of a family 
or Euler-MacLaurin formulae that can be obtained from repeated applications of the 
same argument and, essentially, integration by parts, all to appear in [CS 4]. Our 
earlier approaches used our study of the topology of toric varieties (see [CS 2][S]). 



2. Lattice Points in the Circle 



Let P{t) denote the number of integral points (m, n) in the disk in the plane of 
radius i.e. the number of lattice points satisfying 

+ 71^ < t. 

(2.1)THEOREM. For any e > 0, 

It is known [L] that it is not true that P{t) — ivt is 0{t^) . 

In this section we will derive a new formula for the number of lattice points in a 
circle. 

Let Q > 1 be an odd integer, and let A = q be the polytope spanned by the 
integral points satisfying 

2 / n \ 2 



(9) -'-' 



Let R = ^ and let 



1 ^ ^ 271'i 
fi^^y) = -Q^Yl 1] exp— (px + gy), (1) 

^ p=-R q=-R ^ 



Then 



p{t)= J2fi^^y)- (2) 

Anz2 

In this paper R will be the greatest integer in ^ but the results of this section are 

valid for general Q with the caveat that in (2.2) there would be an additional error of 

t 

It is well known that a circle of radius at most \/t that has any lattice points on it 
has less than 0{t'^) lattice points. Furthermore, for n a positive integer. 
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and it is easy to see that the boundary of A hes between the circle of radius \/i — -q 
and the circle of radius \/t . (For example, if (a, b) is a lattice point in the interior of 
the first quadrant and at most Q\/t — 1 units from the origin, then the triangle with 
vertices (o, 6), (o + 1, 6) (a, & + 1) and the three unit squares with (a, h) as a vertex and 
not containing this triangle all lie within QA .) Hence 

|P(t)-T(A,/)|«^, (3) 

the right side being an estimate on the number of lattice points within a radial distance 
^ of the circle of radius (In the sequel we use the notation f <^ g and |/| < 0{g) 
interchangeably. ) 

A non- horizontal edge Ej has the description Ej — {{Ljy,y)\cj < y < dj} , Lj a 
linear function, and a non-vertical edge the description Ej = {(x, Mjx)\aj < y < bj} . 
Then 

the sum on the right being over the non-horizontal edges. Here and in the sequel 
e(t) = exp(27rit) . Similarly 

over the non-vertical edges, and 



/ ^n5(/-(V'l0^l)-(V'l0 7r2)) - 

^ 1 A 27riq TTin^- f^^ (PLjy + qy. f fr \^ 

E Q2 ^^_^-Q-- .4 <—Q—)My)MLiy) dy 
Ei E ^j''e{'^^^]ML,yKi-5miy 

j ^ P,q = -R 1^ -^"J ^ 

+ Ei E ^re(^-^±J^)^.(M,.)(l-H.))d. 

7 ^ P,q=-R 1^ ^1 -^"^ ^ 
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the first and third sum also over the non-horizontal edges and the second and fourth 
over the non-vertical ones. 
We assert that 



t^+^ InQ 



Q 



We can assume without loss of generality that every edge lies in a quadrant, and 
consider for example the second sum in the first quadrant. Let E be the union of the 
first quadrant edges and M : [a, 6] [c, d] the union of the paramerizations Mj of 
these edges. Then we have to consider, for example, 

EAE E <'^)MMn} = j:^f: Y: <lMn)UQMn), 

3 Q=-flneQZn[a,6] j 5=-^ neZn[a,6] 

with a = ^ , etc., M{x) = Q-^M{Qx) , and ipi{Mn) = ifjiiMn) if (n, Mn) is on the 

interior of an edge and •0i(Mn) = aMj — clmj+i at the indicated vertex, except for the 
obvious modification at a and at b. 

Let N be an integer < N < R + 1 , recall Q = 2R + 1 . Let Pjv consist of those 
n E Z n [d, 6] such there is a nearest integer m to M{n) (there might be two equidistant 
integers) with 

N . r ^1 -/V + 1 

— < I — V rrt -|- n^l < — — — . 

Clearly, 

R+l 

Zn[a,S] = U 



As already noted concerning the number of lattice points in an annular strip of thick- 
ness ^ , the number of elements \Pn\ is 0{{y/t + l)^+^^(5~^) = 0(t2+^(5~^) • Vertical 
distance to a circle is greater than radial distance, so if n G Pn and m is the corre- 
sponding nearest integer to M(n) , then |M(n) — m| > ^^-^ ; this also uses the fact 
noted above that the boundary of A lies within an annulus around the circle of radius 
\/t of thickness ^ . By summation of geometric series, it follows that for n e with 
AT > 2, 
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and this expression is always less than 1. Since \ilJi{QMn)\ < 1, the term we are 
estimating will be at most 



N+l 

\Po\ + \Pi\ + Yl 



\Pn\ t^+'lnQ 



^ N-1 g ' 



which is our assertion. 



Let A = At^Q — Q be the polytope spanned by the points of 

Z[l/Q]eZ[l/Q]n{a;2 + y2 <t} . 

Then, with Ljy = Q^^Lj{Qy) and MjX — Q~^Mj{Q)x , as above, 

T(A, /) = V [ / e{px + qy) + ^ / eipx + qy)MQx) 

+ ~7r" / e(pa; + qy)ipi{Qy) - ^^f^ / e{px + qy)i)i{Qx)i)i{Qy) 
W Ja ^ Ja 

3 ~ 
~J2^T^r7r\ c e{pLjy + qy)iPi{QLjy)dy 

-J^T^T^r-^ / e(px + qMjx)ii;i{QMjx) dx 

~J2^T^^^\ c e{pLjy + qy)iPi{Qy)iPi{QLjy) dy 

- Y^y^P^i a ""{P^ + qMjx)MQx)MQMjx) dx] + 0{^—) . 

3 Q 

(4). 

As Q will less than a positive fractional exponent of t , the factor In Q has been absorbed 
in the exponent of the error. 

Let F{t, Q) be given by the same formula, but with the domain of integration re- 
placed by the disk Df of radius y/t and the boundary with the circle of this radius; 
thus 



12 



SYLVAIN E. CAPPELL AND JULIUS L. SHANESON 



F{t, Q)=y2\ [ eipx + qy) + ^ / e{px + qy)i^^{Qx) 
p,q ^Jd, Q Jd, 

+ j e{px + qy)ipi{Qy) - J e{px + qy)'ipi{Qx)'ipi{Qy) 



Vt 



/ e(pV t-y^ + qy)ipi{QV t - y^)dy 



+ ^/ 4-pVi--y^ + 'iy)i^ii-QVt--y^)dy 

— / e{px + q\/T^-x^)ipi{Q\^T^-x^)d: 

+ — / e{px — q\/ 1 — x'^)ipi(—Q\/ 1 — x'^)d. 
2<5 J-^t 



mq 



Vi 

Vt 

' ■ ^' ' ^' \x 

Vt ^ 

Vt ^ ^ _^ 

Vt 

e{pV t-y'^ + qy)i^i{Qy)i^i{QV t - y'^)dy 



(5) 



VI tiU I 

+ 

e{px + gV t — x'^)'^i{Qx)'^i{Q\/ 1 — x'^)dx 



e{-pVt-y^ + Qy)i^i{Qy)i^i{-QV't - y^)dy 

-Vt 



mp 



+ / e(pa; - q\/t — x'^)ipi{Qx)'ipi{—Qy/t — x'^)dx 

Q J-Vt 

(2.2)PROPOSITION. For any given e > and Q < Vi , 

\T{AJ)-F{t,Q) I « -^In^g. 

Proof of (2.2). We start with the area integrals. For these, let A be the region 
between the boundary ofA on the inside and the boundary of Dt on the outside. 
For j = 0,1 (this actually holds in general) 



Q 



1 ^ 

— I p' exp{2'Kipx) |<^ min{Q, ^} 

p=-R 



where ||a;|| is the distance form x to the nearest integer. As is well-known, this holds 
for j = by summation of a geometric series. For j = 1, 

N N-1 N-j 

p exp{27ripx) = exp{27rijx) exp{27ripx) , 
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from which the result follows using the previous case. It will be enough for the area 
integrals to show that 

^min{Q, ||a;||-i}min{Q, \\y\\-^} < (t^ + t^/^+-^ In^ Q . (*) 

Note that for M an integer and M <a <b < M + 1, 

rb rM+l/Q ^M+l/2 

/ min{Q,\\x\\~^}dx <2 Qdx + 2 — :<ClnQ. 

Ja Jm Jm+1/Q X — M 

Hence, for B a unit square with lattice point vertices, the estimate 

/ min{Q, ||a;||~^} min{(5, lll/ll"^} < / minjQ, min{(5, <C In^ Q • 

JAnB Jb 

Let P be the collection of all unit squares having a vertex within distance ^ of the 
region A. By the same considerations as for (2) above, the number of lattice points 
within this distance of A is at most 

for small (e.g. less than |) e. Therefore, if 



Ic, 



Co^AnljB 

then 

min{g, ||a;||"^} min{g, ||y|r^} < ( ) In^ Q 

'Co \ Q / 

On the other hand, let B a unit square with lattice point vertices, such that AflB ^ 
stays at least a distance from the vertices, a> Q . 
Then we assert that 

/ min{(5, min{(5, lll/ll"^} < ^ ln(5 , (*) 

Let {mB^riB) be the coordinates of the vertex of B closest to the origin. We will 
show why this holds for < ms < ub ■, i-e. the diagonal of B lies in the part of the 
first quadrant above the line y = x , and leave other cases to the reader. 

Let 

Bi = [niB, tob + 1] x[nB + ^,nB + l- 

The vertical thickness of An B over any point on the x-axis is at most since A 

is obviously contained in an annulus of thickness ^ . (Actually the vertical thickness is 
^ ; see below). Therefore, 
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/ min{Q, ||a;||~^} min{(5, |||/||~^} < — / mm{Q, ||a;||~^}(ia; -C — InQ . 

JadBi Q JruB ^ 



li Ar\B C AnBi , the assertion is estabhshed for this square B . If not, then there are 
two possibihties: The region A, moving from left to right, enters this square on the left 
vertical edge, crossing through {uib} x [nB + ^,nB + l — i] and exits through the lower 
horizontal edge, crossing through [ms + ^,mB + 1 — ^] x {ub} or it enters through 
[rriB + ^, niB + 1 - ^] x {ns + 1} and exits through {ms + 1} x [ns + ^, tt-b + 1 - ^] . 
In either case, it is not difficult to see that > 1 • (Look at the positive x- value on 



the circle for y — — {y/t — .) 

Consider the first possibility. It is an excercise to show that the horizontal length 
of ^ n 5 opposite a point on the j/-axis is at most ^ . ( Show that the tangent to the 

circle at a point {z, x) with x > 1 crosses through the line y = z — ^ at a point at most 

a horizontal distance at most ^ from the circle.) Let Ay = An (R x {y} . Let 

B2 = [mB, "^B + ^] X [ub, "-b + ^] • 

The region A cannot intersect [mB,m_B + ^] x [ubjUb + ^] nor can it intersect 
[ruB + 1 - ^] X ["-B, "-B + ^] • Hence 



L 



m.m{Q,\\x\\ ""^l min{(5, ^} < 



AnB2 



nB + - 



min{(5, ||y|| ^} / adxdy <^ ^ In. Q . 

Ja^ q 



Ub 



The second possibility is similar, and this proves the inequality (*) . 

Let A' be the annulus of (radial) thickness ^ with the same outer boundary as A . 
A G A' . For a > b > Q , the number of lattice points with distance from either the inner 
edge or the outer edge of A' in the interval [a~^, b~^] is at most 0{t^^'^'^^{b~^ — a~^)) . 
Hence, if P{a,b) is the collection of unit squares with lattice point vertices with the 
closest vertex to A in this set of lattice points, and 

C(a, b)^An [j B , 

I3{a,b) 



then 



min{(5, ||a;|| ^} min{(5, ||y|| ^}<lnQ 



IC{a,h) 

Let N = N{e) be the largest integer such that Qt~^'^ > 1 . Then 

JV 

^ = Co U U C(Qt-('^-i)% Qt-''') U C{Qt-^\ 1) . 
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Hence 

r ^l/2+2e 

/ min{Q, min{g, \\y\\~^] < N{e) \v? Q . 

J A Q 

With e replaced by 2e , this completes the proof of (2.2) for the area integrals. 
Next we consider the boundary integrals. Let 

^ ^ ^ 27ri ^ n 

\n\ = l 

be the indicated truncation of the Fourier series of ip . Then it is well known that 

Also, the length of dA is 0{^/t) ; this can also be proven using the above remark on 
tangents. This means that we can replace ilJi{±Q^t — y^) with t/j'^^^^Q t — y^) and 
M^QVt^^) with i;^^\±QVt — x^) , in the integrals we are considering, with an 
error no greater that ^ , provided N > Q"^ . (This follows by estimating the exponential 
times iiJi{±Qx) or il:i{±Qy) by unity; by instead summing up a geometric series we 
could reduce A?" to Q at the cost of at this point multplying the error by InQ .) 

Let Ej be a non-horizontal edge. Ej will necessarily lie entirely on one side of the 
y-axis, say on the positive side. Considering one of the boundary terms, for example, 
we have (including the case n = ) 



d 



Q 



e{{p + nQ)^/t - y^ + qy)ipi{Qy)dy - I e{{p + nQ)Ljy + qy)'ipi{Qy)dy 

2'Ki{p + nQ)e({p + nQ)x + qy)il^i {Qy)dxdy . 



Therefore, 



e{pVt^^ + qy)ipi{Qy/t^^)ipi{Qy)dy 

/c (^{P^jy + qy)'ipi{QLjy)'ipi{Qy)dy 



1 /" /t 

V - / 27Ti{p + nQ)e{{p + nQ)x + qy)iJji{Qy) + 0{—) 

In| = l •^^+ ^ 



Here the sum is over edges to the right of the y-axis, c is the greatest rational number 
that is at most \/t and has denominator dividing Q , and is the part of A to the 
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right of the y-axis and with —c<y<c. Notice that the error in replacing \/t with c 
is trivially at most 0(1) . 

The argument for the desired estimate is now the same as for the area integrals 
above, modified as follows: 



R 



R 



R 



^ e{{p+nQ)x) = J2 ^^^^^^ ^ ^e(pa;)+ ^ e{nQx) ^ e{px) 

\n\ = lp=-R \n\ = l ^ p=-R^ 



n =1 



p=-R 



For the first term, 



E 



e{nQx) 



R, 



^e{px) < min{(5, ^jlnQ 



\n\ = l p=-R ^ 



as above except for the extra factor of In Q 
For the second term, within a unit inters 
variation O(QlnQ) . Hence, by integration by parts 



For the second term, within a unit interval we have that "%^/f^Zl[Z'^/l^^ has total 



/ ^ e{nQx) ^ e{px) < ^ — ^ < In^ Q , 

which is what is needed for the argument above for the integrals over Ad Bi . On the 
other hand, on an interval of length ^ the total variation will be at most four times the 
maximum modulus (the derivative changes sign at most twice) and so in the argument 
for ^ n ^2 we have 

„ R 

/ Yl e(nQx) ^ e{px) \ < ^ ^ < ^InQ. 

-^"^y \n\ = l p=-R \n\=l ^ 

Thus we get the same result up to an extra factor of In Q . 

Other boundary terms are dealt with in the same way. This completes the argument 
for (2.2) 



The Fourier series 



, , , 1 sin27ri/t 1 eiut) — ei—vt) 



1 ^(^ 

27ri ^ V 



e{yt) 



can be substituted for ■^i in (5). Then the explicit integration formula 

( 



e{p\/ t-y'^ + qy)dy - I e{-p^/t- y'^ + qy)dy = 27rip / e{px + qy))dxdy 

-Vi J-Vi JDt 
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and the same thing with respect to y can be used to replace the boundary terms by 
area integrals. Some terms cancel, so that what results is the formula 



Fit,Q)=J2 [ <Px + iy)-20 (- + -)/ eiip + fiQ)x+iq + uQ)y) 

P,q=-R "^^t ^ \u\,M=l ^ ^ -^^^ 



+ 2 5Z / ^((^^ + + + ^ / e{px + {q + i^Q)y) 



(6) 



Obviously the third and fourth term are equal when summed over p and q , so that 

n CXD „ 

F{t,Q) = y2\ e{px + qy) + 3y2 e{{p + iiQ)x + qy) 

E {l + l)leiip + ,Q)x+iq + .Q)y)]. 



(7) 



ImI>I=i 



So up to an error of 0{t^'^^Q ^) , this give an expression for the number of lattice 
points in the circle. 

(2.3)Proposition. 

^(*)= E [/ <PX + qy) + zY, / e{{p + iiQ)x + {q + vQ)y) 

p,q=-R "^^t 1/^1=1^* 



(8) 



^ E (- + ^)/ e((p + /.g)x + (g + z.g)y)J+0(tH^g-i). 

I I I I 1 ^ ^ *^ D+ 



A finite expansion can be obtained by substituing ip^^') in (5) for ipi , instead of the 
full Fourier expansion. To estimate the resulting error, we have, as a typical example, 
with (j) either ij) or ij)^^^ , 



eipx + qy)4>{Qy){i;{Qx) - i;^^\Qx)) ^ E 4>{y)e{qy)e{{p + i^Q)x)) 



Vt 



For |p| < i? , by explicit integration 



Hy)e{<iy) 



.e{{p + iJ,Q)x)) dxdy . 



e{{p + iJ,Q)x))dx 



\n\Q 
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It follows that 




e{px + qy)cl>{Qy){il;{Qx)-il;^''\Qx)) « 



Vi 



NQ ■ 



Therefore, with similar considerations with the roles of x and y reversed, the error in 
the fourth term on the right of (5) by using the truncated Fourier series will be at most 
0{^^) . The other terms are dealt with similarly. In particular we have: 

(2.4)Proposition. Let N > . Then 



This is not the best possible result, N — Q is good enough, but an analogue of (2.4) 
will suffice in the sequel. We will not derive the better result in this setting, but will 
illustrate by translating what we have done so far into classical terms. However, the 
rest of this section will play no role in the sequel. 

We recall a little from the theory of Bessel functions. 




(8') 




N 




/ r e(r\/ + 6^ sin(^ + T)d9dr , 
Jo Jo 



with sin r = , ,„ and cos r = , i* ,„ . I 
of Jo (e.g. [W,2.2,(2)]), and [W, 2.12, (5)], 



Hence, by [W, 2.2 (p 20), (7)], the evenness 




e{ax + by) — 2ti 





Thus 
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p,q=-R V ^ y 



oo 



— 1 



(9) 

The term in (9) corresponding to p — q — is the area, 

lim W- Ji(27r\/ta) = Trt. 
a^Q \ a 

In the sum over /j, inside the brackets, we can restrict to < Q , (or < i?) without 
increasing the overall error beyond 0{ti) . For the last sum in the brackets, restricting 
to /i, z/ < Q adds an error of 0(t4 Int). To prove the second statement, for example, 
recall | <^ -It follows that the truncation error will be at worst 

1 „i 'sr^ /I 1\ 1 
t^R^ > - + - 3-. 

To estimate the sum, it sufl&ces to consider the case v > ji . 



H or v>Q 



1 1 1 1 1 InQ 

v>ix>Q tx<Q,v>Q V -w V -w 

which implies the desired estimate. 

Prom just the above, plus the well known expression of Bessel functions in terms 
of trigonometric functions, classical estimates on exponential sums lead quickly to the 
corresponding classical results on the circle problem. For example, the trivial estimates 
yield upon setting Q = t^^^ the result 

In general for {k,l) an exponential pair (see [GK]), we get 

P(t) = 7rt + 0(t^+'). 
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If {k, I, k, I) is an exponential quadruple (conjecturely all pairs give such) then a two- 
dimensional dyadic subdivision argument gives the improvement 

P(t)= 771 + 0(1^^+'). 

For example, the known pair (1/30,26/30) leads by the two methods to the known 
exponents 27/82 + e and 12/37 + e respectively, in the error term. However, in the 
sequel we do not use deep results on exponential pairs. 



3. A Summation, Near the Circle 



We will begin the actual proof of the main result by considering a certain expression. 
Let J denote the set of integers (m, n) with 




^2<n<2Vi. 
The expression is 



m = 



/ / r~2i;(^)(sin^) 2 cos^cos27r(r(\/t csc^ + mcot ^ — n) — |) d^dr . 

(10) 

Here v is a, smooth bump function that is unity for j < 9 < ^ and zero for < ^ < ^ 



and ^ < 9 < TT . Most results to follow will be true for general Q and R but when 

needed it will be assumed R = [t^"*"^] • 

We will expand the inner integral / = /(r, m, n) using the stationary phase result 
as given in [H, 7.7.5]. The simplest way to do this is to set 

CSC 9 = cosh X cot 9 = sinh x. 

Then d9 = sechx dx and 

/oo ^ g 

vi{x)r~'^ tanha;sech2 (x) cos 27r{r{^/t — m? cosh (a; — A) — n) — ^)dx . 
-oo 

In this expression, vi{x) = v{9) and 

A \/i + m 



Vt — w? 
Clearly 
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r~^vi{x + A) tanh(a; + ^) sech^ (a; + ^) cos 2^T{r{^/T--rr^? cosh x — n) — ^)dx . 

-oo 

From [H, 7.7.5], / will be the real part / = 3?/i of 
h = ^^''^^^ ~ T ""^^ V r-^{t - rn^yi V r^2-'^L»2.-3M(tanh-(secli)i)(-yl) 

j = 0(2) 2i^>3ji 



+ 0(r-^+^(t-m2)-^). 

(11) 

(There are only even j in the summation because odd derivatives of cosh vanish at 
zero.) In particular, 



/= 6 cos2T^r{^t-m^ - n) + 0{r~^ {t - rrT)) . (12) 



As a consequence 

;2 rR 



S{t)= V Xll^ f cos27Tr{^/r^ - n)— + 0{\nH) . (13) 
We will need some summation results related to the expression that appears in (13). 



(3.1) Proposition. Let ry > . Let a be the greatest integer in \/tj2 . Let L = L(t, 77) 
denote be the region in the plane with < x < a and 



\y 



- Vt-X^\ < =T 



Then 

Vt — m? ^ , r. n N dr 

(m,n)€L 



Evt — f „ , / . ( 
Y / cos27rr{-\/t — m'^ — n)- 

4^ r sin2T,{rlt^l)%+t--^C{R,t) + 0{t'^+'). 



Notes: 1. The result is stated and proved for R = 1^4+"^] . In general, the error 



would have an extra factor of . 

2. The essential feature of the term C{R,t) is its independence of r] . However, we 
will actually obtain the following expression for this term [N > or even N > Q , 
see (2.4)): 

C{t,R)= Y: E r^"^27r((p + A^Q).. + gVF^)^^^^ 
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The proof of (3.1) will involve the application of the following form of the classical 
Euler-MacLaurin formula: 

J2 0(P)= [\iz)dz+ [\{zyi;{z)dz + i:{a)cj>{a)-mm ■ 

a<p<b 

Here 

The first integral on the right will be referred to as the principal term, the second as 
the remainder term, and last two expressions as the endpoint terms. 
Let L' differ from L by removal of the curve y — \Jt — — r . Let 

dr 



E= y. / cos{2nrF{m,n))- 



r 



(m,n)£L' 

with 

Fix,y) = y- \Jt-x^ . 

Then it is not difficult to see that t~3E differs from the sum of (2.5) that we are 
studying by at most 0{t'^) . 

More generally, let -{t"'} - I < di < d2 < n + 1 and let L{di,d2) be the set of 
points {x, y) in the plane satisfying the above condition on x and 



dr 



di < y — \/t — x'^ < d2 ■ 

Let 

S((ii,(i2)= / cos(27rrF(m, n)) 

{m,n)EL{di,d2) ^ 

We now apply the above form of the classical Euler-MacLaurin formula to y- 
summation in this expression. The endpoint terms obtained in this way and summed 
over m provide a contribution to 'E{di,d2) of 



^i^iy) J cos{2nrF{m,y)) 



dr 
r 



m 

This will be considered later. 

The principal term in j/-summation produces a contribution 



nVt-m'^+di r-R 

y2 , / cos{27rrF{m,y))—dy . 



It is easy to sec that the endpoint terms for this x- summation will amount to no more 
than 0(\n.R{d2 — di)) , and the principal term will be 
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^ dr 
cos{27TrF{x,y)) — dxdy . 



This will also be considered later, and we will now look at the remainder term for this 
x-summation, which will be (with a change of order of integration) 

27r / / ip{x)Fx{x,y) sm{27rrF{x,y))dxdydr . 

Jl JL{di,d2) 

The part of L{di, ^2) not contained in the cone given in polar coordinates by 



tan-i^^ = ^o<^<- 

is contained in a triangle of area at most {d2 — cZi)^ . Since as well as the other 
factors of the inner integral have modulus at most one, it follows that removing this 
region affects the entire expression by an error of at most R{di — ^2)^ • On the other 
hand, the boundary curves of the intersection of L{di, (^2) with this cone have the polar 
coordinate equations 

/ d'^ 

p = disinO + Jt - dlcos9 = Vi + di sin6' + O(^) . 

Therefore, again as the modulus of the integrand is less than unity, if K is described 
in polar co-ordinates by 



and 



tan-^^^ = ^o<^<J 
a 2 



\/t + di sin 6 < p < \ft + 6^2 sin Q 



then replacement of L{d\^ (^2) by X affects the integral by at most an error a constant 
muptiple of -R(max(if) . Therefore, we may consider instead 

/"f /•Vt+d2sin6' 19 COS ^ 

III •0(pcos6') e{p sin 6 — \/t — p"^ cos"^ 9) pdpdOdr , 

switching to complex coordinates. 

However, in the domain of integration. 



3 



. ^ / n ^ ^ / cot^^csc^ , ^,9 ^/maxldii ^ 

psin^ - y/t - p2 cos2 9 = csc9 ■ {p- Vi) + ^ (p - Vif +0{ ^) 



and 
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pcos$ . _/max|(iik 
^ cote + 0{ ^) 



Therefore, up to errors of 

0[R^{d2- di)m.3x.{\di\^,d^]t~^) and 0(i?(d2 - c^i) max |cZj|) , 
it suffices to consider instead 



sin 6 

/ / / 'ip{p COS 9) cot 9e{r CSC 9- {p-Vi)+rt~^ cot'^ 9 CSC 9- {p-Vt)'^)pdpd9dr. 

sin 



Up to an error of at most a constant multiple of R{d2 — di) (using p y/i and 
9o > -j), this is the same as 



R r§ rd2sine 



Vi / 'i/'((v^ + 2;) cos6')(cot 6')e(r2; esc 6' + rt ^ coi^ 9 esc 9)dzd9 dr. 

Jl JOo Jdi sine 

To estimate this integral, we plug in the Fourier expansion 

1 eipx) - e(-px) 



2m ^ p 



Therefore, we can look at 

^ rR rd2Sine 



/ / / (cot6')e(r2;csc6' + rt~^2;^cot^6'csc6' + »(2; + v^)cos6')d^ci6'dr. 

l^^27ripA ^ 4 sin/ 

However, the 6'-derivative of the phase function is large (note f < ^ < f in the range 
of integration); in fact. 



I - CSC 6* cot 6* - rt 2 2;^ (2 cot 6* csc^ 9 + cot^6'csc6') +p(\/t + 2;) sin6'| > \p\^t . 

Therefore, by integration by parts with respect to 9 , followed by trivial estimates with 
respect to the other variables, this sum will be at most a constant multiple of 



00 ^ 

R{d2 - cJi) ^ - 



p2 



Hence the remainder term for y-summation of the principal terms for x-summation is 
estimated as at most a constant multiple of 
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R{d2-di)+R{d2-di)'^+R'^{d2-di)maix{\di\^,d'^}t ^ +R{d2-di)meix\di\+Rmeixd'^ 



For e and r] small enough, e.g., e + 2r] < | , and \di\ <S , this will be at worst 
Next, since Fy = 1 , the remainder term from y-summation will be 



/'\/t—m'-^+d2 

— 27r^^ / / 'il>{y) sm{2nrF{m,y))dydr 

= — 271 / / •(/'(a/ t — + z) sm{2'Krz)dzdr 



m 

Summing up, so far we have 



/"^ dr 
S((ii,(i2) = 2_,'^(^) / cos(27rrF(m,2/)) — 

r dr 
+ / cos(27rrF(a;, J/)) — dxdy 

+ 27r^ ^MH^!±^(cos27ri?^ - cos27r^)d2 + 0{t-^+'+^\d2 - d^)) 



= ZTT > / (cos zttRz — COS 2nz)dz . 



Vt—m'^+di 
Vt—rri^+d2 



m 



Note that the first term on the right vanishes if di = —d2 ■ 

To analyze the last term in the above expression for S((ii, ^2) , the remainder term 
from y-summation, let 



A{z) = / {\/t - + z)dx 
Jo 



be the area of the region A{z) under the curve y = \/t — x^ + z and let L{z) be the 
number of lattice points in the region bounded by this curve, the x-axis, and the lines 
X = and x = a . Then, by trapezoidal approximation, 

A{z) = +J2Wt-m'' + z) + + 0(1) , 

m=l 



and 



L{z) = J2 - ™^ + z} + a + l. 



m=0 
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Hence 



J2 ^iVt - m? + 2) = A{z) - L{z) + ^{a + Vi+ \/t-a^ + 2z) + Oil) 



We now indicate the application of a slight variant of our method of sections 1 and 
2 to study the number of lattice points L{z) in the closed region A{z) bounded by the 
curve 



y 



= yjt — + z , 



the lines x = , x = a and the x-axis. The notation of section 1 will be used. Let 
w he a weighting of the lattice points in A{z) given as the limit as the radius goes to 
zero of the proportion of area of a small circle around the lattice point that lies within 
A{z) . Let 

Let 

R 



p,q=-R 



Q 



Then 



L'{Z)^ [ f-{do7Ti) -{60712). 
J A 



As in the proof of (1.1), 



/ /•(5e07ri)-ao7r2)= / /-/ /-(l-^eOTTi)- / / • (1 - 5^ O TTa) 

Jqa{z) Jqa(z) Jqa{z) Jqa{z) 

+ [ /■(l-5e07ri)-(l-5,07r2). 

jQA(z) 



By integration by parts, 



/ / ■ (1 - OTTi) = - / fx-'ipl,e°T^l + 

Jqa{z) Jqa{z) Jo 



Q{Vt^+z) 



-QiVt+z) 

+ I f{x,y)lpl,e{x) 
Q{Vt^+z) 



f{x,y)i;i,e{x) 

^QH-{y-Qzy 



Qa 



dy 



x=0 



dy. 



x=0 
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The second integral vanishes and the third one evaluated at the lower limit is 0{lnQ) 
Therefore 



Jqa{z) „ o ^ Ja 



p,q=-R 



e{px + qy)i(;i^^{Qx) 



Q JA{z) 

rVi+z 

+ 7^ / e(pVt -{y- zf + qy) i^iAQV^ - iv - ^Y)dy 



+ 0{lnQ). 



Similarly 



p,q = -R 



e{px + qy)'4^i,e{Qy) 



Q Ja{z) 

+ — / e{px + q{\/t -x'^ + z) ipi^eiQiVt - a;2 + z))dx 
Q Jo 

Multiplying the ^-integration by parts that leads to this with 1 — Se o tti , we have 



/ / ■ (1 - 5e OTTi) • (1 - 5e 7r2) = - / /y • V'l.e °7r2 • (1 - 5e OTTi) 

jQAiz) JQA{z) 

+ r f{x, y/t-x''+ z)iPiAQ^/QH -x^ + Qz){l - 6,{x))dx . 
Jo 



Similarly to the above, integrating the first term on the right by parts with respect to 

X, 



I /•(l-5e07ri)-(l-5e07r2) = 

JQAiz^ 

An^pq 



R 

E 

p,q=-R 

27Tiq 



A{z) 



e{px + qy)'ipi,e{Qx)'ipi,e{Qy) 



rVt+z 

+ / e{px + g(\/ 1- x'^ + z))il^iAQ{\/ 1 - x"^ + z)dx 
Q Jo 



- / f{x, - ^2 + z) i^l,e{\/QH - X2 + Qz) 5e{x) dx + O(lnt) . 

Jo 

Similarly to the polytope case, upon taking the limit as e ^ 0+ , the expression 
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r f{x, Vt-x^ + z) MVQ^t -x^ + Qz) 5,{x) dx 
Jo 

will become 

E eiqiVr^ + z) MQiVt^ + z) . 

^ m=0 q=-R 

Here ipi{Q{Vt - + z) = V'i(<5(VF — + z) for < m < a , and at the endpoints 
of the range of summation the value will still be at most | in absolute value. 

We will now estimate the contribution of this term to the remainder term from 
y-summation that is under study, i.e., we want to estimate 



— / e{q{\/t — w? + z) ipi{Q{\/t — rri'^ + z){cos27tRz — cos27r2;) — . 

Q ^ n D^di ^ 



m=0 Q=-fl"^'^i 

By summation of geometric series, 

R 



J2 <(l{Vt -m? + z) < mm{Q,\\^/t - m? + z)\\-^} , 

q=-R 



and 

cos 27r Rz — cos 27rz 



<min{g, 1^1 

z 

Hence our contribution will be at most a constant multiple of 

E=^y2 [ ' min{Q, \\Vt - + ^)|r^}min{Q, \z\-^}dz . 

For —R — 1 < N < R , lei Pn denote that set of integers me [0,0;] for which there 
exists an integer m with 

N ^ r ^ + 1 

— < V ™ + — < — — • 

Each m is contained in least one such set; hence 

^-AEE / niin{Q, II - m2 + min{Q, \z\-^}dz . 

The vertical distance of a lattice point to a circle is larger than the difference of the 
distance of the point to the origin and the radius of the circle, and for any interval of 
integration of unit length or less, it follows similarly to considerations in section 2 that 
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/ min{|^ — A| ^, Q} min{|;2| ^TQ}dz< 

J a 



A 



From this and the fact that the number of elements of Pn is 0{t^'^^Q ^) ■, it follows 
that for iV > 2 , 



f ' min{Q, II - m2 + min{Q, \z\-'^}dz < \N - l\-H^+'+'^ 



On the other hand, for jATj < 2 , 



^ / min{Q, ||Vt-m2 + 2;)||-^}min{Q, |2;|-^}d2;< — — / min{Q^ 2;-^} ^2; 



Q vy_^ Jo z'^j 



Therefore, absorbing a factor of In Q into , 



We have 



L{z) = L'{z) + + a + Vt - + 2z) + 0(1) 



and it is also not hard to see that errors 0(1) contribute at most O(lnQ) to the remain- 
der terms from j/-summation under study. Letting e — > O""" in the above expressions, 
we then get 
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I {cos 2t:Rz — cos2t:z) dz — — I I 2. I e{px + qy) 

(P, 9)5^(0,0) 

R 

+ E 

p,q=-R 

27rig 



2nip 

Q Ja{z) 



e{px + qy)i(;i{Qx) 



+ 
+ 



Q JA{z) 
An'^pq r 

Ja{z) 



e{px + qy)ipi{Qy) 

e{px + qy)ipi{Qx)ipi{Qy) 



/ e{p^t -{y- zY + qy)i^x{Ql{^t-{y-zY))i,x{Qly) dy 

J s/t-oi^^-z 



\dz 



As in section 2, ipi can be replaced by the partial Fourier expansion 

2m 



1 



without increasing the error. By explicit integrations, 



rVi+z 

/ e{{p + iiQ)^fr^{y^^ + qy)dy 



p-\- nQ 



P + /iQ 



I e{{p + iJ,Q)x + q{\/T^^x^ + z) dx 
Jo 

ret pVt+z 

/ e{{p + ijlQ)x + g(\/ 1 - ^- z) dx + I e{qy) dy , 

Jo JVt^+z 



for g . The last two terms are 0{q ^) , by integration by parts. Similarly, 



rVi+z 

/ e{{p + iiQ)^t -{y- zf + (g + vQ)y) dy 



P + fiQ 
+ 



^ f e{{p + ijQ)x + (g + iyQ){Vt - x'' + z) dx + 0{{q + vQ)'^] 
Jo 
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and, for any p, q with q , 

[ e(px + qy) = / e(px + 1 — x'^ + z)) dx + 0(q~^) 

Ja{z) 27rzg Jo 

Finally, also for any p , 

fda r rd2 



n f , ^cos27tRz-cos27tz , r / X , f \ ^ „ ^ s , 

/ / e{px) az = / e{px) ax / [coszttHz — coszttz) az 

Jdi Ja{z) ^ Jo Jdi 

r , X I ^ , f"" COS l-nRz - cos 27r^ , 

+ / e[px)\/t — x^dx / az . 

Jo Jdi ^ 



The first term on the right always vanishes, as does the second, under the assump- 
tion di = —d2 ■ Hence, upon subsituting the partial Fourier series and making this 
assumption, we obtain 



> / (cos zttRz — COS 2nz)dz = 

rd2 f ^ '^^1 

\ J2 J2 - e{{p + fxQ)x + q{Vt^+z))dx 



- E E Y. OyiXov\ \ + i^Q)^ + (9 + yQ)^\/t - a;2 + z)) dx 



COS ^'kRz — COS 21^ z , ^ / 1 2 +^+^ ^ , . 
dz^O{ +QlnQ) 

2^ 



It now follows that for d^ = —di and di = 7 e Z , 



y I (cos 271 it2; — cos zTTZJa^ 

^ Jdi ^ 



m 



1 fh+'^+n 

= — i3(7) + o(— — + ging) 



with 
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B(7) 



/ cos 27r((p + + q\/ 1 — x^)dx ■ / 

Q 



sin 27rg2; 



52; 



cos 2'kRz — cos 27rz 



jr- / cos27r((p + + + yQ)\/t - x^)dx 

p,q 1*^1=0 
kl=i 



'sin27r(g + z/(5)-2 



{q + i^Q)z 



cos 27ri22; — cos 2ttz 



dz . 



In general, 



sin27rA2; 



cos 27r uzdz 
z ^2 



1 n sin27r;2(A + /x)) , 1 F sin27r^(A - /i)) 



-dz 



;(Si(27r(A + /.)7) + Si(27r(A-/.)7) 



Hence, 



with 



^(7)= 2 (^+(7) + ^-(7)), 



R 



B±(7)= E E 



' COS 27r((p + iJ,Q)x + qyjt — x"^) 



dx- 



Si(27r(Q ± i2)7) - Si(27r(g ± 1) 



R 



^ ^ Qiy Jo 

p,q=-R lfl=o ^ ^ 



g r COS 2tx{{p + + + pQ)\/t - x^) 



g + i^Q 

Si(27r(g + vQ ± R)-f) - Si(27r(g + z^g ± 1) 



dx- 



The function Si(x) is odd, bounded, and, for positive x , 

Si(a;) = ^ + 0(max{l,a;-^}). 

It therefore foUows that 



SOME PROBLEMS IN NUMBER THEORY I: THE CIRCLE PROBLEM 



33 



with 



P,q7^0 |m|=0 



B(7) = C{t, R) + ^(7) 



q2 /J- 

/"V 2 COS 2tt{{p + nQ)x + q\/t - x^) 



dx, 



and 



R 



' COS 27r((p + /i(5)a; + g\/t — x'^) 



dx 



+ 



+ 



^\q\ 7max{l, — Qf} 'ym.ax.{l,R + q} 



R 

p,q=-R 1/^1=0 ^ 
1/1 = 1 



' COS 27r((p + //(5)x + + i'Q)Vt - x"^) 
q-\-yQ 



dx 



1 1 

+ 



7|z/Q + Q + i?| 7|z/Q + Q-i?| 
It is clear that i?) is indeed independent of r] . 

For general p > q> 1 , the derivative of the phase function px + qy/t — x^ satisfies 



X 



27r(p- ?^==)| :^p-q. 



Vt-x'^' 



Therefore, 



/ COS 27T {px + q \/t — x'^) d: 
Jo 

Similary, if p < g < — 1 , 

/ COS 27r {px + q \/ 1 — x"^) dx 
Jo 

whereas if p and q have opposite signs, then 

/ COS 27r {px + q \J t — x^) dx 
Jo 



p-q 



1 



Finally, the absolute value of the second derivative of the phase function will be 



\q\t 
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and hence 



i-a. 

I COS 27r ipx + q \/ 1 — x"^) dx 
Jo 



t4 



These estimates imply |i3(7)| <^ 7 -"^t 4 In i?, as we shall now explain. 
For example, 



1 

— J cos 27r{px + q^/ 1 — x'^)dx 

— / cos 27t{px + q\J t — x'^)dx 

R>\p\>\q\>^^ 

+ — / COS 2'k{px + q\/ 1 — x'^)dx 



p,q = — R 



U 



« E 



<Clni? + t4 < 



Similarly, but using only the second derivative estimate. 



R 

E 

p,q = — R 
q^O 



\q\ max{i? ± g, 1} 



pa 

/ COS 27t{px + q\/ 1 — x'^)da 
Jo 



R 



« E 



< 



+ 



p,«=-« ig|*niax{i?±g,l} |g|^i? V kl niax{i? ± |g|, 1}. 



This takes care of the terms in the first summation in the above estimate for B{'y) with 
11 = . For the remaining terms, I < \fJ,\ < ■, 



R N 



\^ I f ^ cos27r {{p + iiQ)x + qVt^^) 



R N 



1 



q^O 

R N 



•c 



E E l^|Qg2+(^Q)2(|^| + |p + ^Q|_|^|)«l^^- 
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Also, 



E E 



p,q = -R 



^ IqI maxji? ± Q, 1} 
m|=i ^ ^' 



dx 



R N 

« E E 



q^O 

R N 



max{i?±g,l}|/x|g(|g| + |p + ^g|-|g|) 



p,q= 

q^O 



For the other summation, for the terms with fu. = , we can use only the second 
estimate involving the second derivative to get an estimate of 0{t^Q~^) . This is left 
to the reader. Fo the remaining terms, we use the first derivative estimate on the 
integral when |//| > l^-] and the second derivative estimate for < l^-] . This leads to 
an estimate, up to a constant, of 

R N ^ 



E E 



p.^-M ,,t;:',.,M\<i + ''Q\(.\p+i'Q\-\t + ''Q\KWQ + t\-R) 

ImI>I''I 

R N 1 N 1, ^ 

+ V V T- ^ < In^ Q + V — ^ < 1 4 

p,q=-RM.M=i W\\q + i^Q\HWQ + q\-R) 1^1=1 |/i|^<5^ 

ImI<I-I 



This completes the estimation of ^(7) . 

What we have so far can be summed up as follows, with 7 = r = {f^J : 

E{-T,T) = jC{t,R) + 2 f f cos27Tr{F{x,y))—dxdy + 0{t^+^'^+'). 

The integral can be replaced by one over the same region as in the discussion of the 
remainder term for x-summation above. Because of the factor of ^ in the integrand, 
the quadratic term in the Taylor series of the phase function is not needed, and we 
obtain the following integral: 

fR rVt+T sine 

/ / / cos 2n {r {Vt — p) CSC 6) pdpdd — 

Jl J 60 J Vi-T sine 

fR I.+T sine , 

= Vi / / cos27r{rzcsc9)dzd9— + 0{Rt'^) 

Jl Jeo J-TsinO 



\ft dr 
j- / / (sin27r(rr) - (sin27r(-rr))sin^d^^ + 0(i?t'^^ 



a f 

Jl 



1 Jeo 

^ dr 
sm27rrT — + OiRt"^) 
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This completes the proof of (3.1). 

There is the expansion, due to integration by parts, with aU non-zero coefficients, 

dr ^ ^ 



•^1 n=0 



sin27r«;r^ = «2n+i«;-'"+' + 0{U-^) + 0(«;-2^-i) 



Using this and the preceeding, we may restate the previous result as follows (ab- 
sorbing some constants into the Oj): 



(3.2) Proposition. Let r] > Let a be the greatest integer in Yt/2 . Let L — L{t, a) 
denote the be the region in the plane with < x < a and 



Then 



/ cos27rr( Vt - -n)— =t-^C(t,R) 



i=0 



To summarize so far, let 

J = Ji U J2, , 

with (m, n) G jTi if (m, n) & J' and 

\Vt-m^-n\ <T= in 
and (m, n) e J^2 otherwise. Let 

S{t) =Si{t)+S2{t) + 0{lnH) 

be the corresponding decomposition of (13). 
Then, from (13) and (3.2) 

Er{t) = t-"^C{t, R) + t-^J2 «2i+iT-''-' + 0(t''+2^) . (14) 

i=0 

4. Away From the Circle 

To study £2 1 split J2 into J2 and J2 , consisting of those (m, n) e J2 inside and 

outside the circle of radius , respectively. We will use complex notation, specifically 
replacing cos27r(-) with e(-) . for the term being considered. By integration by parts 
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72 rR 



N 



+ 7V! 



\/t — 



1 (27rV^)^n^r^+i(Vt^^-n)^+ 



(r\/ 1 — m^) 
-e{r{\/ 1 — w? — n))dr 



(15) 



For the term involving the integral on the right, the trivial estimate gives 



R 



Vt 



1 1 n^r^+^{Vt^^ -n)^ 
which will be at most 0{t~^/'^~^n~^) for 

AT = N{r]) > 



e(r(-\/ 1 — m? — n))dr 



< 



Int 



n 



4ri 



Hence the sum over of this term will be at worst 0(1) • 

We can rewrite the sum over J'^ and over J'{' of the i*'* term as 



ce / no(m) 



Vt — 'W? 



(27rV^)V^+i^^o 



n 



{Vt — w? — ny 



I e{r\/ 1 — w?) 



R 



and 



Vt — m? 



(^-1)! y-/ 



I e(r\/ 1 — m^) 



respectively, where a = 1^/^/2] , r?,o(m) = |Vt — m^ — r] , and ni(m) = |Vt — m^+r] . 

If the inner sums were monotone functions of m , Abel summation and Lemma 3.5 
of [GK] or even just Kusmin-Landau would apply as above to estimate these terms. 
This type of argument does work if the inner sums are replaced by integrals. Namely, 
consider 



(^-1)! 
(277 V^Yr'+^Vi 



VT- 



Vt — w? du 



Vt 

2 



u 



{Vt — m? — uY 



I e{r-V t — m?) 



and 



(z-1)! / /-Pv^J 

(27r V^) V^+i Vt^X Jvt^+r 



Vt - m2 du 
u{u — Vt — m?) 



- I e{—r\/ 1 — m?) 



R 



38 



SYLVAIN E. CAPPELL AND JULIUS L. SHANESON 



Then it is not hard to show that the integrals are monotone functions of m , just 
differentiate them with respect to m. By Kusmin-Landau or [GK,3.5], 

a 

1 

and, more generally, 

1 

It follows by Abel summation that the first expression will in modulus be at most a 
constant times 

Jl u{Vi-uy Ji \U y/i-uj 

and similarly, the second will be at most O(lnt) as well. 

To study the error introduced by replacing sums by integrals, we use a one di- 
mensional Euler-Maclaurin expansion, which we now review. There is a unique con- 
tinuous periodic funtion tl)2{t) of unit period, with V'2(^) — V'i(^) fo'^ u ^ Z and 
with ifj2{u)du = 0. Then define il)q{u) inductively as the unique continious func- 
tion satisying ip'j{u) = ifjj-iiu) and Jq iljq{u)du = 0. Given real numbers a < 6, let 
X(w) = X[a,6](w) = 1 for u e (a, b) , x(w) = for u ^ [a, b] , and x(o) = x{b) = | . Then 

Yx{n)<P{n) = / <P{u)du + Y,i-iy{Ub)<l>^'~'Hb)-i;,{a)<l>^^-'\a)) 
nez "^^ j=l 

l-b 

J a 

For example, we apply this to the first sum. Note that 



d^- 



Vt — m? 



du^ ^ u{\/t - m? - uY 



z(z + l)---(z+j-2) 



u{Vt — — uy+^' 



In particular, when evaluated at u = \\lt ^ this is 0(t 2 Also, for N = Ni = 
N{ri)-i, 



d 



N 



Vt 



du^ u{Vt -m? - uY 



Finally, there are at most O(t^) values of m , if any, such that \/t — m? is also an 
integer, i.e. there are at most this number of lattice points on the circle of radius \/t . 
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It follows that the error in replacing the sum with the integral in the i*^ term, for the 
region J'+ , will be 



[i-iy. 



+ Oif). 

Since 



Vt — m?'^j ( \/t — m?) 



elryt — m^) 



y/t — w? — T y/t — m? 



for < m < q; , without increasing the error we can replace this with 



(i-1)! 



a Ni 



t . . . (i + - 2)^i<^H^e(.y^) i; + 0(t') 



m=Q j=l 



Similarly, the error for , will be 



(-1)' 



V V / V ffi=0 T — 1 



a N, 



m=0 j 



However, it is clear that for r = R, these terms are 0{t^R ^) , so finally we can 
consider just the sum of 



a Ni 



and 



(-1)^ 



[i-iy 



We have (compare [R]) the Fourier expansions 



and 



(-1)'= ^ e{pu)+e{-pu) 
^^^-^^"^ = (2.)--^ ,f ^ 



(271) 



lpl=l 



p 



2k 
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The series are absolutely convergent from -02 on, and by Kusmin-Landau or [GK, 3.5] 
for j > 1 odd our error contributions will be constant multiples of 



+1^1 gt^v^TT 



and will actually vanish completely if i and j have different parity. For j = 2k > 2 
even and i also even , however, we will have a contribution to the error; we get, for a 
given i , 



{27rV^y ^ (27r) V^+^-i ^ ^ ^ " 



j=2k 



For j = 1 and z even the terms cancel. For j = 1 and i odd, we get a contribution of 



? y 



m=0 

a 



(2,rv/=T)'+'T'(J 5. 



E 

|p|=l 

oo 

E 

lpl=l 



e((l - TO^) 



P 



e{\Jt — m^) — 1 



e{pVt — rn? ) e{p\/t — m? 

^ ^ p{p + 1) ^ ^ 



p=2 



P(P - 1) 



The argument using Kusmin-Landau or [GK] then appplies again to the appropriate 
terms, so finally, for a given i , we end up with 



r + O(r-0 • 



We give what is probably a well known interpretation of the last sum. By trapezoidal 
approximation using integer points from to ck , 



where A is the area under the circle of radius y = yt — x^ for < a; < 6 and L the 
number of lattice points in this regions, counting the top and sides but not the bottom. 
Combining this with "Pick's theorem" for the number of lattice points in the triangle 
with vertices the origin, (q;,0) , and (a, a) , it follows that 



^^(V^) = ^ - L. + ^ - + 0(l) , 
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Li the number of lattice points contained in the sector of the circle of radius with 
angle f < ^ < f . However, it is obvious that 

P(t) = 8Li + 4y^-4Vt, 

and hence 



8 

i 

To sum up, the have shown, with a slight change of notation that the part £{t) 
corresponding to the portion of J that is far from the circle has the form 

JV(77) Ar(,7) 

E2{t) =t-^Yl ^2i+lT-'^-' + t-^TTt - P{t)) J2 C2i+lT-^'-'0{f) . (16) 
i=0 i=0 

Explicit summation expressions can be written down for the various coefficients, and 
it is clear that in particular ci ^ . 

5. Conclusion of the argument 

The combination of (14) and (17) gives 



N{rj) 

N{r)) 

i=0 



(14) 



Recall that in this equation r = T{r]) = . Now subtract this equation from the 
same result for ^ We can take N{^) — N{r]) or else let the extra coefficients be zero. 
Also, multipy hy . This yields 



i=l 

Since ci 7^ , the sum in the right side is invertible for large enough t and 
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It is known that irt — P{t) changes sign for arbitrarily large t ([L, VIII, 6 or 7]; see also 
[T, 12.6] and [H, 1916a]). Therefore, the coefficients a2i+i + d2i+i must vanish at least 
until the terms not of higher order than the remainder; this vanishing could also be 
verified directly. In any case. 



follows. 
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